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, Q ■ ABSTRACT 

^— i , A simple model for the evolution of disc galaxies is presented. We adopt three numbers 

■ from observations of the Milky Way disc, Ed the local surface mass density, the 

stellar scale length (of the assumedly exponential disc) and v c the amplitude of the 

£f) . (assumedly flat) rotation curve, and physically, the (local) dynamical Kennicutt star 

' formation prescription, standard chemical evolution equations assuming a Salpeter 

IMF and a model for spectral evolution of stellar populations. We can determine the 
0^ ■ detailed evolution of the model with only the addition of standard cosmological scalings 

with time of the dimensional parameters. A surprising wealth of detailed specifications 
follows from this prescription including the gaseous infall rate as a function of radius 
and time, the distribution of stellar ages and metallicities with time and radius, surface 
brightness profiles at different wavelengths, colours etc. Some of the detailed properties 
are as follows: the global gas infall rate and the global star formation rate are almost 
constant at 2—2>M ( 7 ) yr~ 1 and 2-4M |/r~ 1 during the evolution of the disc. The present 
O-f day total mass in stars and in gas is 2.7xl0 10 M© and 9.5 x 10 9 M Q , respectively, and the 

q • disc has an absolute K-band magnitude of —23.2. The present day stellar scale length 

?-h ' (normalised to 3kpc) in the K-band and is larger than at shorter wavelengths. At the 

, solar neighbourhood stars start to form w lOGyrs ago at an increasing rate peaking 4 

C$ ■ billion years ago and then slowly declining in good agreement with observations. The 

mean age of long lived stars at the solar neighbourhood is about AGyrs. The local 
• i-H , surface density of the stars and gas are 35 and 15M Q pc -2 , respectively. The metallicity 

■ distribution of the stars at the solar radius is narrow with a peak at [Z/Zq] = —0.1. 

The present day metallicity gradient is — 0.046dexfcpc _1 and has been significantly 
steeper in the past. Using a Chabrier IMF increases the luminosity of the model 
and results in a steeper metallicity gradient. The local metallicity distribution is only 
weakly affetcted. Different formulations for threshold densities for star formation have 
been tested and lead to a truncation of the stellar disc at about Ylkpc. Comparisons 
with the current and local fossil evidence provides support for the model which can 
then be used to assess other local disc galaxies, the evolution of disc galaxies in deep 
optical surveys and also for theoretical investigations such as simulations of merging 
disc galaxies. 

Key words: methods: analytical - galaxies:disc - galaxies: formation - galaxies: 
evolution - galaxies: fundamental parameters 
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1 INTRODUCTION 

At the present time our knowledge of the current proper- 
ties of nearby galaxies, including, especially our own Milky 
Way, far exceeds our ability to model the prior evolution 
of these systems. Working forwards in time starting with a 
reasonable cosmological model and computing via ab initio 
methods the evolution at high spatial resolution is a goal 
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being pursued by many rese arch teams (e.g. lAbadi et alJ 
l2003al : iRobertson et al.ll2004l and references therein) and it 
may ultimately be successful. However, the project is hand- 
icapped by two potentially crippling inadequacies: insuffi- 
cient numerical resolution and inadequate representation of 
important physical processes. We shall attempt in this paper 
to follow an alternative approach which is essentially the in- 
verse of the paradigm described above. We will take as given 
the present gross properties of our Galaxy: the disc surface 
mass distribution is exponential and the rotation curve is 
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flat. Using these and other properties of normal disc galax- 
ies, we will see if we can construct a prior evolution that 
ends in this state and is consistent with our knowledge of 
cosmology. Then we can compare the results with our de- 
tailed knowledge of the Galaxy, the distribution of stellar 
ages and metallicities, to see if our hypothetical evolution- 
ary path is consistent with the present observed state. 

The dark matter, largely assembled into self gravitating 
halos, we take as primary cosmological fact. The distribution 
of halo masses roughly follows the Press- Schechter distribu- 
tion 
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where n is the sp ectral index of the power spectrum 
jLacev fc Coldll994D . 

Thus, to a good approximation, the evolution of the 
halo distribution can be computed when the evolution of the 
nonlinear mass scale m„ (t) is known. This in turn is fixed by 
the evolution of the linear power spectrum of perturbations 
Pn{i) oc K n since 
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defines the nonlinear wave number knl and m»(£) = 
const. < p(t) > /itm at any give n n. 

Following lLacev fc Cole) (Il994f) we note that halos much 
more massive (at some point in time) than m,(t) are rare 
and are being assembled by unusual merger events in iso- 
lated high density regions and their number density is 
rapidly increasing with time. But halos having m « m,(t) 
have been left behind. They live in lower density regions (or 
else they would have been merged out of existence) and their 
numbers are slowly declining. 

At the present time the nonlinear halo mass scale is of 
the order of 1O 14 M0 but individual halos within which or- 
dinary galaxies live are of order 1O 12 M0. Thus, these have 
either been merged into bigger systems (i.e. these halos are 
now hosts of 'cluster galaxies') which tend to be bulge domi- 
nated, or else they live in a lower density 'field' environment 
and have prominent disc components. Hence if we consider 
the halo of a system such as our own Milky Way the most 
likely history was for it to have been assembled via mergers 
until a time tf when its mass reached « 1O 12 M corre- 
sponding to m*(£f). Then it was, for whatever reason (most 
likely the less than typical density of the larger region sur- 
rounding it) , side tracked from the cosmic merger dance and 
survived, with a modest level of subsequent accretion, to the 
present day. Thus, we take it to have been embedded in a 
halo with virial velocity v v [ T which increased with the cosmic 
typical values v v i T (t) = v* il: (t) during the merger phase of its 
evolution until reaching its current value at tf after which 
the galaxy entered its quiet evolutionary phase. The virial 
velocity remained essentially constant at its present value. 
The halo, however, still has grown in size and therefore in 
mass. During the quiet evolutionary phase the gas within 
the halo can cool and add to the central disc component 
thereby increasing the inner circular velocity v c (r,t) due to 
the increasing disc mass. For simplicity we take v c (r, t) to be 



flat, independent of radius, in accordance with observations 
of typical spirals and the idea that barygenic infall increases 
the mass density in the inner parts of galaxies to approxi- 
mately _th^l^vel^Hh£maximum of the dark matter rotation 
curve llFlores et al.lll993l) . 

Since disc galaxies now have close to exponential pro- 
files and there e xist several reasonable arguments for how 
this came about llMestell963l : lFall fc Efstathioull98d: iGunnl 



1982; Dal canton et al.Nl997l: ISlvz et alJ 120021) . we will as- 
sume that the total surface mass distribution was always 
exponential with a scale length r<i(i) scaling with the cir- 
cular velocity as given by spherical infall theory r^(t) oc 
Vvir(t) / H(t). This simple prescription suffices to specify the 
infall of gas to the disc at all times in the past as a func- 
tion of both time and radius. Then, if we couple this infall 
rate with a model for s t ar for mation abstracted from ob- 
servations by iKennicutd <ll998T) . that at each radius 10% of 
the gas is turned to stars per rotation period, and add a 
standard treatment of chemical evolution, we have a com- 
plete prescription for the evolution of the stellar and gaseous 
components at every radius. It is amazing that this simple 
model is able to successfully fit a large number of observa- 
tional data. 

Thus, to summarise, our assumptions are that, after 
an initial formation period, the disc always maintains an 
exponential profile with fixed central surface density in a 
halo providing a flat rotation velocity and the star formation 
follows a Kennicutt-law. The three numbers known from the 
present, rd(to),T,o(to),v c (to) fix the final state with v v ir(£) 
and £p(£) held const ant and r^(t) scaling with 1/H(t) (see 
iFerguson et alJ 120041 for observational evidence). At early 
epochs (t < tt), when the scaling relation would give an un- 
typically big and massive system (for the epoch) , we simply 
follow scaling laws of the typical halo at these earlier times. 

Over the last decade chemical evolution models have 
been used to explore the formation of the Galaxy and 
successfully reproduce properties of the solar neighbour- 
hood, abundance gradients, metallici tity distributions, 
photo me tric properties etc. (e. e [ Beckman fc Pa gel! 
19891: iPardi fc Ferrinil Il994l : iPagel fc Tautyaisienel 



199a IPrantzos fc Aubertl Il99f ; IChiappini et alJ Il997l: 
Molla etalJ Il997t lAllen et alJ Il998t IPrantzos fc Silkl 
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Chang et alJ Il999t iHou et alJ \200A iRomano et al 
Ma^uc£^fe"Cbd^3plnj 
Chang et alJ 120021: ICasuso fc Beckmanl 12004 . Many of 




these models define the infall rate as a function of time 
and radius, e.g. exponential in time with e-folding time 
scale s increasing with radius (see e.g. IBoissier fc Prantzos 
^99[) or Gaussian-shaped infall rates IpTantaos^^Biik 
1998) to fit the local properties like gas surface density, 
abundance gradients or the metallicity distribution. Some 
authors have used cosmological scaling relations to asses 
the possible range of present day disc properties. However, 
the cosmological evolution of the discs with time, which is 
governed by the infall rate (mostly exponential), has been 
fixed to fit the present day properties. The infall rates were 
not directly connected to the cosmological evolu tion (see 
e.g. Ijimenez et al]|l998t IBoissier fc PrantzodEoOoll . 

In this paper we follow a different path. We assume 
that the evolution of the infall rate onto the galaxy is fully 
constrained by the cosmological model and the small set 
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of boundary conditions described above and investigate the 
properties of the resulting galaxy. Those properties are then 
compared to the Milky Way. We are aware that such a model 
that is aimed to reproduce global properties of early type 
discs will probably not be able to reproduce all detailed 
observations our Galaxy like e.g. colour that depends on 
the very recent star formation history. On the other side 
as we consider the Milky Way being a typical spiral galaxy 
the model has to be able to reproduce its global properties. 
Given the small number of initial parameters and physical 
processes involved it is surprising how well the model can 
reproduce Milky Way properties. 

The paper is structured as follows: A review of mea- 
surement of the properties of the Milky Way that are used 
to scale our model are given in Section|!2] In Section|3]we de- 
scribe the details of the model. We develop the scaling rela- 
tions that determine the infall rate, describe the star forma- 
tion rule used and the way we have implemented chemical 
evolution. Our model predictions on the assembly history, 
the metallicity, and photometric properties are compared to 
the Milky Way properties in Section^] The effect of changing 
the IMF form Salpeter to Chabrier is investigated in Section 
|S| In Section HJ we investigate the influence of a threshold 
value for star formation on the outer profiles of the disc and 
in Section we discuss our results and conclude. 



2 PROPERTIES OF THE MILKY WAY 

The model for disc galaxy evolution which will be described 
in Section |21 needs three input parameters: the present day 
scale length of the total disc surface density distribution 
rd(to), the central surface density of the disc Ed(io), and 
the circular velocity v c (to) (assumed to be flat). In this pa- 
per we restrict ourselves to a comparison with the Milky 
Way. This approach has both advantages and disadvantages. 
If we consider the Milky Way as a typical early type spiral 
galaxy, any successful model should be able to reproduce the 
gross properties of its disc. In particular the model has to be 
able to reproduce data from the solar neighbourhood that 
mainly depend on the integrated past star formation like the 
gas and stellar surface mass density and the metallicities of 
the gas and the stars. Other properties like the colour will 
directly depend on the very recent star formation history 
and might not be typical for the mean properties of early 
type spiral galaxies. Measurements of global properties of 
the Milky Way like total disc mass and total luminosity suf- 
fer from large observational errors. Here the model can help 
to constrain the data given that it is in agreement with the 
solar neighbourhood observations. 

In this Section we will review the latest observational 
constraints on the three parameters that we will later use 
as the input for our model: the scale length of the stellar 
disc rd,*, the circular velocity at the solar neighbourhood v c 
and the total surface density of the disc Ed (to) at the solar 
radius r . 

2.1 Distribution of stars and gas 

The stellar disc of the Milky Way is assumed to follow an 
exponential surface density profile over a wide range in ra- 
dius. Measurements of the scale length are in the range 



rd,» ~ 2.5 — 3.5fcpc (see Sackettl jl997f) for a review). More 
recent investigations report values above as well as below 
3fcpc, e.g. star counts from the 2MASS s urvey yield ~ 3.3fcpc 
llLopez-Corredoira et alj|2002h whereas IZheng et all (1200 ll) 
get ~ 2.75kpc based on direct HST observations of M stars. 
For the model presented in this paper we therefore assume 
a conservative value of rd,» = 3fcpc. Several authors find 
the exponential profile does not fit the outer parts of the 
Galaxy. They find an edg e (truncated stellar emissivity) 
at 4 — 6fcpc from the Su n jRobin et al.lll992t iRuphv et alJ 
ll996UFreudenreichlll99gl) . However^ lLope^Corredoln^t^dl 
(2002) find no sign for a cut-off in the stellar disc at radii 
r < 15fcpc. Observations from other disc galaxies clearly 
indicate an abrupt change in the slope of the stellar pro- 
file in more than 60% of the observed disc g alaxies at radii 
of about 3 — 4 times the disc scale le ngth ([van der Kruitl 
[ 19791: Ivan der Kruit fc Searld Il981allpl: iPohlen et alJ l2000l : 
Ide Griis et aUl200lft . It has also bee n proposed that t he 
galactic disc might have an inner hole dRobin et al.ll2003h . 

For the dista nce of the Sun from the galactic centre we 
adopt r — 8kpc llReidll993D and for the loca l circular speed 
we use I)© = 210fcms _1 (see Sackctt (f997) for a review). 
Both values are slightly lower but still in agreement with 
the IAU standards for the galactocentric distance of the Sun 
and the local circular velocity of rp, = 8.5 ± 1. lfcpc and 
v Q = 220 ± 20fcms _1 jKerr fc Lvnden-Beljll986T) . 

The radial distribution of atomic and molecular gas in 
the M ilky Way is shown in Fig. The data from iDamd 
( 1993) have been kindly provided by Thomas Dame. The H2 
is more confined to the centre than the HI which can extend 
to large radii at almost constant surface densities. This is a 
global feature that is observ ed for external disc galaxies as 
well llYoung fc Scovirldll99lD . The central depression in HI 
seems to be common for early-type type disc galaxies like 
the Milky Way or Andromeda, whereas only about half of 
them shows a central depression in H2 llYoung fc Scovilld 
Il99ll) . We also note that the Gal actic ./^-distribution shows 
a prominent peak at the centre (ISanders et al]|l984[) . 

2.2 Local surface density 

The most reliable measurements have been made for the to- 
tal mass density w i thin about lfcpc of the Galactic plane. 
iKuiiken fc Gilmord dl989t Il99ll) used spatial and kinemat- 
ical data of local dwarf stars to determine the vertical 
gravitational potential and the total mass density (disc 
+halo) to be 71±6M pc -2 within 1. lfcpc from the Gala ctic 
plane near the Sun. Recently iHolmberg fc Flvnnl (120041) re- 
ported a similar value for the total dynamical mass within 
1. lfcpc of 74±6M pc~ 2 , 6 5±6M<t)PC~ 2 within 0.8 kpc, and 
41M pc~ 2 within 0.35 kpc. lSiebert et alJ (I2003T) determined 
the total mass density wit hin 0.8fcpc to be 76l 2 2M pc -2 
and lKorchagin et alJ «2003f) determined the dynamical mass 
within 0.35fcpc to be 42 ± 6M pc" 2 . 

The determination of which fraction of the matter 
resides in a flattened disc component or in a flattened 
or spheroidal halo depends on further assumptions, e.g 
the disc scale length, the scale height, the galactocentric 
dist ance of the Sun and the distribution of dark mat- 
ter. IKuiiken fc Gilmord Jl99lf) found a local disc surface 
mass density of 48 ± 9M pc -2 from modelling the Galac- 
tic rotation curve. Some values determined by other stud- 
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ies depending on ad ditional modelling are 841. ^ Mq pc 



feahcall et al.1 ll992Tl . 52 ± 13M F) pc- 2 iFlynn fc Fuchsl 
1994h. or (Y7 + *lMvpc- 2 jSiebert et alJ l2003h . Recently. 
Holmberg fc Flvnnl l20olh found 56±6M pc" 2 for their disc 
model. 

In contrast to dynamical measurements of the disc sur- 
face density, direct observations of the stars at the solar 
neighbour hood should be less dep endent on addition al as- 
sumptions. |Couideiai] dl996l) and lZheng et al.l fcOOll) used 
HST observations of M stars in the Galactic disc to di- 
rectly determine a surface density of 12.2 — 14.3Mqpc~ 2 . 
Visible stars other than M stars contribute « 15Mqpc~ 2 , 
resulting in a total stellar surface density in the range of 
~ 27 — 30A/qpc -2 . Alternatively, star counts within 5pc 
of th e Sun using Hipparcos parallaxes djahreifi fc Wielenl 
I1997D combined w ith the observed vertical disc profile 
( Zhenfl et al.||200l h giv e a similar number of 26.9Mqpc~ 2 
ijBinnev fc Evans! 1200 IT) . With a local gas surface density 
(the one we assume h ere is slightly higher th an the one de- 
rived bvlDamelll993l) of » 13 - 14M pc~ 2 jBahcall et alJ 
Il992l: lOlling fc MerrifieldlEooih the total mass density of 
the disc would lie in the range of 40 — 44M pc _2 . Inter- 
estingly, those measurements give surface densities that lie 
significantly below the dynamical est imates. As the effe ct of 
binary stars is supposed to be small (iGould et al]ll996l) the 
difference might be caused by undetected dead stellar rem- 
nants and brown dwarfs. These we estimate to 5 — 7M pc -2 , 
to give a total surface density of 45 — 51M@pc~ 2 . 

Independent of which number is closer to the truth there 
is a discrepancy between the amount of identified disc mat- 
ter and the total column density within l.lkpc of about 
20 — 35Af pc~ 2 which could be attributed to a spheroidal 
dark matter component. The exact amount of dark mat- 
ter in the solar neighbourhood is unknown. Most authors 
find no evidence for dark matter in the disc at the sola r 
radius jKuiiken fc GilmoreilT99lt iHolmberg fc Flvnnll2004l) . 
(However, taking the dynamical measurements and the di- 
rect observations at face value 10 — 15M pc~ 2 could reside 
in a dark disc component, see next Section) 

Throughout the paper our model is constrained to have 
a total mass surface density of the baryonic disc components 
(stars + gas) within \z\ < l.lkpc at the solar radius r = 
8kpc of 50Mqpc~ 2 . The division between live stars, stellar 
remnants and gas is determined by the computation. 



where the central surface density £0 and the scale length r d 
change with time. The cumulative mass distribution of the 
disc is 

M d (r,t) = M d (t)(l - (1 + r/r d ) exp(-r/r d )), (4) 
where 



M d (t) = 27r£ (f)rd(f) 



(5) 



is the total mass of the disc. In the limit of a thin disc its 
rotational velocity is given by 

v 2 d (r,t) =4^GSoWr d (t)j/(r,t) 2 [Io(!/)ifo(y)-/i(y)ifi(y)],(6) 

where G is the gravitational constant, y(r,t) = r/(2r d (i)) 
and Ii(y) and Ki(y) are the m odified Besscl functions of the 
first and second kind jFreemanll970h . Real discs have a finite 
thickness and a slightly lower rotation velocity. The effect is, 
however, small and will not affect the results presented here. 
The rotation velocity of the disc peaks at r^.iit) = 2.15r d (i) 
at a value of 



< 2 .2 it) = 0.774ttGEo (t)r d (t) 



(7) 



with y = 1.075 in Eqn. HJ The disc material within T2.2 
accounts to 64.5% of the total mass of the disc (Eqn.[IJ. At 
7"2.2(£) the disc contributes a fraction of 



/v(t) =V d ,2. 2 (t)/v c (t) 



(8) 



to the circular velocity v c (t) which is assumed to be con- 
stant with radius. In general, f v can vary with time. The 
exact value of f v for present day disc galaxies is still under 
debate. If galactic discs were maximum then f v — 0.85 ±0.1 
at r2.2 which s eems to be close to correct for the Milky Way 
( Sackett 1997). However, u sing data from nearby disc galax- 
ies lCourteau fc Rlxl il999l) favour a value of f v — 0.6 ±0.1 
which is in good agreem ent with the value of f v = 0.63 de- 
rived bv lBottemal 1^993). Our adopted model has a present 
day value of f v = 0.61 (see Fig.Hjl. 

Using the above scalings we can estimate the total mass 
in dark matter (assuming a spherical distribution, e.g. the 
dark matter does not concentrate significantly in the disc) 
inside T2.2 using 
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(9) 



3 THE MODEL FOR DISC GALAXY 
EVOLUTION 

In this Section we describe the theoretical framework we 
need to model the formation and evolution of an early type 
galactic disc. 
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(10) 



which gives 4.3 x 10 10 M© (1.9 x 1O 1O M ) for /„ = 0.6 
(/„ = 0.85). Assuming a constant d ensity of a spherical dark 
matter halo up to the solar radius feinnev fc Evans! <l200ll) 
argue that no halo with a cusp steeper than r~ 0,3 is viable) 
we can estimate the column density within zi.i = l.lkpc 
using 



3.1 Disc scalings 

We assume that in the absence of star formation the gas in a 
given halo would settle in a disc with an exponential surface 
density 



£ d (r,t) = £ (t) exp(-r/r d (t)), 



(3) 



Sd m (r2.: 
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(11) 



With T2.2 = 6-Gkpc, v c = 210fcms _1 and f v = 0.6 this 
would result in a dark matter surface density of 26.4M0pc -2 
or 11.3M pc" 2 for the maximum disc case (/„ = 0.85) 
and leave 45 — 60Mp)PC~ 2 for the disc itself, assuming the 
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Table 1. Predicted local properties of model and observations of the Milky Way at the solar neighbourhood 



Observable 


Model 


Observed values 


Reference 


Total surface density from stellar dynamics 


\z\ < l.lkpc, [M & pc~ 2 ] 




71 ±6 


Kuiiken fc Oilmore (19911 






74 ±6 


Holmberg & Flvnn f 20041 


\z\ < 0.8kpc 




' u -12 


Siebert et al. (20031 






65 ±6 


Holmberg & Flvnn (20041 


\z\ < 0.35fcpc 




42 ±6 


Korchagin et al. ("2003") 






41 


Holmbera- & Flvnn ("20041 


Total disc surface density from disc modelling 


[M @ pc- 2 ] 


50 


46 ±9 


Kuiiken & Oilmore (19891 






48 ±9 


Kuiiken & Gilmore (19911 






84+29 


Bahcall et al. (19921 






52 ± 13 


Flvnn & Fuchs (19941 






56 ±6 


Holmberg & Flvnn (20041 


Surface density of visible stars from direct observations 


[M pc- 2 ] 


32 


35 ± 5 


Gilmore et al. (19891 






27 


Gould et al. (19961 






30 


Zhena et al. (20011 


Surface densities of stellar remnants 


[M @ pc- 2 ] 


3 


2-4 


Mera et al. (19981 


Gas surface densities from direct observations 


[M @ pc- 2 ] 


15 


8 ± 5 


Dame (19931 






13 - 14 


Oiling fc Merrifield (20011 


Star formation rate 


[Mepc-iGyr- 1 ] 


6.4 


2-10 


Guesten fc Mezger (19821 






3.5-5 


Rana (19911 


Infall rate 


[MQpc^Gyr- 1 ] 


3.3 






Luminosities 


L b [L B qPC~ 2 ] 


38 


20 ± 2 


van der Kruit (19861 


L v [L v@ pc- 2 ] 


30 


22.5 ±3 


Pagel (19971 


L k [L K qPc~ 2 ] 


59 


68 ±23 


Kent et al. (19911 



Table 2. Some observed global properties of the Milky Way and model predictions 



Observable 


Model 


Observed values 


Reference 


Star formation rate [MQj/r ] 


3.6 


0.8 - 13 


Rana (1991) for references 






3.5-5 


" Rana (19911 


Infall rate [A/qj/t -1 ] 


2.2 


0.5 - 5 


Braun & Thilker (20041 


Total luminosities 


Lb[Lb@] 


2.4 


1.8 ±3 x 10 10 


van der Kruit (19861 


l v[Lvq] 


1.9 


(1.4 < L ve < 2.1) x 10 10 


van der Kruit (19861 






1.24 x 10 10 


Bahcall & Soneira (19801 






2.1 x 10 10 


Sackett (19971 


Lk[Lk@] 


4.1 


4.9 x 10 10 


Kent et al. (19911 
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3.2 Cosmological scalings 

During the self-similar structure formation process, the time 
evolution of the halo depends on t he assumed cosmol- 
ogy. In the spher i cal collapse model llGunn fc Gotd Il972t 
lBertschingedll985l : ICole fc Lacevlll996ll the virial radius r v j r 
and the virial mass M v i r (the mean density within r v i r is 
200/Ocrit) of a halo with virial velocity v v ir at any time t is 
given approximately by 



Vyix(t) _ 

10H(t) ' 



M vil = 



10GH{t) 



(12) 



rvir(t) = 



with 

H[z(t)] = J ffo[nA,o + (l-nA,o-no)(l+^) 2 +n (l+^) 3 ] 1/2 , (13) 

where H(z ) is th e Hubble parameter at redshift z (see also 
iMao et al.1 Jl998l) 1. We chose h = 0.70 and fi = 0.30 for a 
close d universe as constrained recently by WMAP and SDSS 
data llTeemark fc et aLlEooll) . 

We assume that the formation and evolution of a typ- 
ical disc galaxy proceeds in two phases. At early cosmic 
times the galaxy follows the general cosmological evolution, 
e.g. we have estimated the evolution of the circular velocity 
by calculating the redshift at which the value of the Press- 
Schechter mass function is maximum for a given halo mass 
(Antonio Vale, private communication). In other words, we 
take the peak circular velocity of the most common halo 
at every epoch. At the time a halo reaches its present day 
virial velocity the second phase starts which is the domi- 
nant phase of disc growth. We assume that the halo decou- 
ples from the general cosmological evolution and its virial 
velocity will stay constant from then on: 

v v ir(t < tform) = Wvir(f) of most common halo (14) 

V v ii(t ^ tform) = ^vir (tform). (15) 

We call this time the formation time tform of the halo. For 
simplicity we assume that the virial velocity of the galaxy is 
flat at all radii. 

We now assume that the scale length rd of the gas disc 
forming within a given halo is a fixed fraction f T of the virial 
radius r v i r of the halo, 



Td(t < tfor 111 } 
m J 



/ r r- v ir(*) = /i 



Vyirjt) 

10H{t) 



(16) 



r ^vir(tform) , ..-^(tform) / n \ 

J r — 77iTr77\ — = r d (tform) — TTT7\ — • (I') 



WH(t) ■»v-»™> H (t) 

After its formation time the virial radius of the halo and 
its total mass do still increase according to Eqn. 1121 The 
circular velocity, v c , at the inner parts of the halo is now 
corrected for the disc component that is building up, 



){t > tform) = «vir(iform) + «d,2.2(*)- 



(18) 



For simplicity we do not take the effect of the adiabatic 
contraction of the dark matter halo into account. 

In the upper panel of Fig. Q we show the evolution of 
the circular velocity of the disc galaxy. Initially it follows 
the evolution of the most common halo at every epoch. At 
the formation time, tf olm « 2.5Gyrs, the halo decouples 
from the cosmological evolution and the circular velocity 
increases due to growth of the disc to its present day value 
of 210fcms~ 1 . As shown in the second panel of Fig.^the disc 
scale length evolves according to Eqn. 1161 and 1171 assuming 



ft = 1/70. This value has been chosen to guarantuee the the 
final stellar scale length is 3kpc. The scaling implies that the 
total disc angular momentum J(t) oc Mdt> c fd scales with the 
expected angular momentum of the infalling gas. 

Combining Eqns. 171 ISl and HrJl we obtain an expression 
for the evolution of the central surface density of the disc. 
At times earlier than tform we keep f v — const, and let 
the central surface density increase; after tform we keep the 
central surface density fixed as the infalling gas will have 
a too large angular momentum to fall to the centre. As a 
result, f v will begin to increase as f v oc H~ 1 ^ 2 : 



E (t< tform) - J 
/v(t < tform) = COnSt., 



10tf(t)Vvir(t), 



and 

Eo(t ^ tfo; 



So (t* 



/rEo(tfi 



(19) 
(20) 

(21) 
(22) 



fAt> M = °- 7 ^ G WH{ t) Vc (ty 

In the third panel of Fig. Q we show the evolution of the 
central surface density of the disc. The system has been 
scaled according to a present day total disc surface density 
of 50MqPC~ 2 at r = 8kpc, resulting in f v = 0.36 at tf or m (see 
Section |2| for the parameters assumed for the Milky Way). 
For times t < tform, fv has been kept fixed at this value, at 
t ^ tform, fv is increasing according to Eqn. 1221 to a present 
day value of fv = 0.61 (see last panel in Fig. [Q. 

By now the total surface density profile of the disc is 
known at any time (Eqn.|3J. Therefore the surface density 
of the infall rate is given by 

E d (r, t + dt) - E d (r,t) 



(r,t) 



fit 



(23) 



To summarise, the evolution of the infall rate of the disc 
is fully constrained by three present day parameters: the 
circular velocity of the disc, the disc scale length (defined 
by / r ), and the central surface density of the disc (defining 

fv at tform)- 



3.3 Star formation law 

At this point we have to include a model for star formation to 
the model to follow the evolution of the stellar and gaseous 
phase separately. After the gas within a ha l o has settled into 
a disc it starts forming stars. iKennicuttl 1^998) has given 
two alternative formulations that successfully describe the 
relation between the averaged surface density of the star 
formation rate (SFR) and the averaged total gas surface 
density in observed disc galaxies. The first can empirically 
be parametrised by a simple Schmidt-type law 



(24) 



where on average n ~ 1.4. The SFR surface density of indi- 
vidual galaxies can, however, deviate by a factor of 7 from 
this relation. The second formulation, which seems, empiri- 
cally, equally valid, and is more physically motivated, is 

Esfr oc ^Si, (25) 

Tdyn 

where Td yn is the dynamical time of the system (oc f2 -1 ). 
A law of this kind was predicted in case the passage of spi- 
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Figure 1. Upper panel: Circular velocity v c of the model galaxy 
versus cosmic time (solid line). The dashed line shows the peak 
circular velocity of the most common halo at every epoch. The 
galaxy halo reaches its final circular velocity 190fcms _1 ) at 
the formation time tf orm ~ 2.5Gyrs (indicated by an arrow in 
the lowest panel). Thereafter v c increases due to the growth of 
the disc to its final value of v c = 210kms~ 1 . Second panel: Evo- 
lution of the scale length with cosmic time according to Eqns. 
Hoi and IT7I assuming f T = 1/70 (solid line). The dashed line for 
comparison shows the evolution for r oc H~ . Third panel: Total 
central surface density of the disc versus time according to Eqn. 
1191 and 1211 scaled to the Milky Way value at t = if orm and kept 
constant thereafter. Lowest panel: Time evolution of / v fEqns. l^ul 
and 1551. 



ral arms plays an i mportant role triggering star formation 
jWvse fc SilkHl989t) . In the original formulation, iKennicuttl 
( 1998) found that on average 21% of the available total gas 
mass per orbit within the outer edge of the disc is trans- 
formed into stars, or 



7~orb 



(26) 



where T or b = 27rr out / v c is the orbital p eriod at the outer ra- 
dius r out and e w 0.1 iKennicuttll998t) is the star formation 
efficiency. Alternatively this relation can be written as 



where r gas = E gas /EsFR is the time scale for the depletion 
of gas by star formation. 

Throughout this paper we use a formulation based on 
the local dynamical time (rotation period) of the system. At 
every radius the surface density of the star formation rate is 
given by 



Esfr(t, t) 
with 



T OI b(r, t) 



Torb 



(r,t) = 



27rr 
v c (t) 



(28) 



(29) 



and a star formation efficiency of e = 0.1. IWong fc Blitzl 
(2002) have tested the validity of the above local formu- 
lations, which were originally derived from disc averaged 
data, for individual discs using azimuthally averaged data 
of a sample of molecule-rich spirals. They find that the 
parametrisations based on a local Schmidt-law and on the lo- 
cal dynamica l time of the systems are both consistent with 
observations. iBoissier et alJ {2003) who i n addition tested 
one f urther parameterisation proposed by iDopita fc Rvdeil 
1994) 



(30) 



where Et is the total surface density, that is sometimes used 
in chemical evolution models C Matteucci fc ChiappinilboOll) 
find a good agreement with the Milky Way data for all three 
formulations (howe ver with different value s for the expo- 
nents m and n than Do pita fc Rvderl (I1994h ). In particular, 
they find a good agreement assuming a flat rotation curve 
for the Milky Way. 

It is not surprising that the above simple prescriptions 
for star formation that depend on the gas density itself fit the 
observed data reasonably well. All algorithms which main- 
tain a relatively moderate ratio of gas to stars will lead to 
the same results. 



3.4 Chemical evolution 

We follow the evolution of the model galaxies in independent 
rings assuming no radial gas flow s using a modified version o f 
the chemical evolution model of lOstriker fc Tinslevl dl975l) . 
In every independent ring the change in gas surface density 
Eg and surface density in stars E s is given by 

dS e {r,t) = -T, SFR {r,t)dt + K- ms {r,t)dt 

+if latc (r, t)dt + Ei FR (r, t)dt (31) 
dE B [r,t) = E S FB.(r,t)dt - K ins {r,t)dt 

-K ute (r,t)dt, (32) 

where Esfr is the star formation rate per unit area (Eqn. 
1281 and Eifr is the rate of gas infall onto the galaxy per 
unit area, as defined in Eqn. 1281 K[ ns is the mass per unit 
are in gas ejected from massive stars instantaneously, Ki&te 
is the mass per unit area in gas ejected at later evolutionary 
phases of low mass stars. They are defined as 



K ins (r,t) = i? ins E SF R(r,t), 



Torb 



(27) 



#late(r,t) = 



t (t',r)W(t-t')dt'. 



(33) 
(34) 



8 Thorsten Naab and Jeremiah P. Ostriker 




Redshift 
4 3 2 1.5 1 0.6 0.3 



0.1 



Time [GyrsJ 



Figure 2. Fraction of mass lost from a single burst population 
of stars versus time according to Eg ns . 15*31 and 1331 for a Salpeter 
IMF (solid line) whic h is a good app r oxima tions to the metal de- 
pendent value of lBruzual fc Chariot! i2003t) shown by the dashed 
lines for 0.0001 < Z < 0.05. 



Table 3. Input parameters for the Milky Way disc model 



Input parameters 

Solar radius 
Circular velocity at tq 
Total disc surface density at tq 
Disc collapse fraction 
Star formation IMF 
Star formation efficiency 
Solar metallicity 

Effective yield 
Z of infalling gas 



tq = 8kpc 
v c = 2Wkms- 1 
S Q = 50M Q pc- 2 
fr = 1/70 
Salpeter 
e = 0.1 
Z Q = 0.0126 
Y = 0.0095 - 0.0135 
= 1 x 10" 4 - 3.7 x 10" 



Rins = 0.1 is the fraction of gas returned instantaneously to 
the ISM from newborn massive stars and W(t) is a weighting 
function defined as 



W(t) = R, 



5, - 1 

TO 



t + T 



(35) 



with R„ = 0.3 , <5* = 1.36, and tq = 1 x 10 s assuming a 
ISalpeterl Il955ft IMF ICiotti et aljil99ll) . This analytic ex- 
pression is a good approximation (see Fig. |5J for the frac- 
tion of returned gas for a single burst population to the 
metal dependent values o f the spectral evolution model by 
iBruzual fc Charlod teOOSt) that we have used to compute the 
photometric properties of the model disc (see Section l4.4|l . 
The evolution of the mass of metals in the gas z g follows 

E g (r,i)dz g = Y(l - z g {r,t))X S FT<.{r,t)dt 



+ (zw 



,(r,t))T, WR (r,t)dt 



,(r, *') - z g (r, t)]Es FR (r, t')W(t - t')dt'dt 



(36) 



where Y, the yield, is the mass fraction of newly formed 
metals of the total mass of gas returned to the ISM. zif is 
the metallicity of the infalling gas (see Section f4.3> . 




Figure 3. Upper panel: Time evolution of the total mass, the 
mass in stars and the mass in gas for the model disc. Lower panel: 
Time evolution of the total surface density, the stellar surface 
density and the gas surface density at the solar radius. 



4 MODEL PROPERTIES AND COMPARISON 
TO THE MILKY WAY 

In this section we present the properties of the galactic 
model disc assuming the present day properties of the Milky 
Way discussed in Section [5] and the evolutionary model of 
Section |3 All input parameters are summarised in Table |2] 
In addition, we compare the model disc to available global 
and local observed properties of the Milky Way. 



4.1 Mass assembly 

The upper panel in Fig. [3] shows the evolution of the total 
disc mass in gas and stars, respectively. The total mass of 
the disc out to 26kpc is M tot = 3.7 x 10 10 M Q with about 
M g = 9.5 x 1O 9 M in gas and M tot = 2.7 x 10 10 M© in 
stars. These values are in g ood agreement with mass mod- 
els for the Milky Way by iDehnen fc Binnevl dl99gl) that 
are based on observational constraints. For r<j = 8kpc 
and rd/rQ — 0.375 they would predict a total disc mass 
(ISM + thick stellar disc + thin stellar disc) in the range 
of 3 — 4.5 x 10 10 M(7) depending on additional theoretical 
constraints. iBahcall fc Soneiral lll98CT) estimated the stellar 
mass of the disc to be 2.0 x 10 10 M Q . 

The time evolution of the local surface density is shown 
in the lower plot in Fig. |3 The model has been normalised 
to a present day total surface density of 50Mqpc -2 . The gas 
starts to assemble at the solar radius after 3Gyrs, and after 
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Figure 4. Stellar surface density profiles (upper panel) and gas 
surface density profiles (lower panel) at 2, 4, 6, 8, 10, 12, and 13.6 
Gyrs (as indicated in Fig.[3]for the solar radius). The present day 
profile is shown by the thick lines. 



GGyrs the local gas surface density stays almost constant at 
its present day value of 15M G pc~ 2 in good agreement with 
the observed 13 — \<±Mqpc~ 2 . The first stars at the solar 
radius form after AGyrs followed by a steady increase to the 
present day value of 35Mqpc~ 2 (with « 3Mqpc~ 2 invisible 
in stellar remnants) resulting in « 32MqPC~ 2 visible stars 
which is only slightly higher than the value of 27— 3QMqpc~ 2 
derived from direct observations of normal stars (see Section 
I2.2l for discussion). 

The time evolution of the surface density profile for 
stars and gas is shown in Fig. 01 and reveals the inside 
out formation of the stellar disc, originally suggested by 
iMatteucci fc Francois! 1^389). The stellar profile is always 
close to exponential with increasing scale lengths with time. 
The gas profile develops a central depression with an expo- 
nential outer profile. The present day profiles are indicated 
by thick lines. In Fig. |S] we show in addition the present 
day total surface density profile, which is an exponential by 
definition. The stellar profile is very close to exponential in 
the outer parts, however, with a smaller scale length. The 
gas is almost depleted at the centre. The surface density 
is rapidly increasing with radius and reaches a maximum 
at Ikpc. Thereafter it decreases almost exponentially with 
a scale length significantly larger than the one of the stars 
and of the total surface density. The gas distribution of the 
model disc compares well with the observed one (lower panel 
in Fig. 01. The peak, however, is at a slightly larger radius 
(Ikpc compared to the observed bkpc). We have to note that 



Figure 5. Upper plot: Final surface density distributions of the 
total baryonic matter (dashed), stars (solid) and gas (dotted). 
Lower plot: Observed surface density distribution for H2 (dotted 
thin), HI (dot-dashed thin), and H2 + HI (solid thin) increased 
by a factor of 1.4 to account for Helium and heavier elements. The 
data has been kindly provided by Thomas Dame <Dameill993T) . 
The thick solid line is the model prediction for the total gas den- 
sity. The thick dotted (thick dot-dashed) line indicat es the model 
predi ction for H2 (HI) assuming Shi/Sjj 2 oc r~ 1 iWong Kr, Blitd 

we did not explicitly include a central bulge (bar) compo- 
nent to the model which might have a strong influence on the 
evolution of the gas profile in the inner parts, e.g. star forma- 
tion close to the centre might be suppressed and the Milky 
Way bar might induce gas inflow towards the centre. There- 
fore the profile at small radii is not to be over-interpreted. It 
is, however, surprising that we get a reasonable fit without 
a detailed modelling of the bulge. Of course, by assuming 
a flat rotation curve into the centre we are, dynamically, 
implicitly allowing for a steep bulge component since the 
rotation due to the disc alone would be linear (vd oc r). 

4.2 Infall and star formation history 

The time evolution of the radial distribution of the surface 
density of the star formation rate and the infall rate (see 
Eqn, l23H for the modelled galactic disc leading to the inside 
out formation is shown in Fig. |S| At early times the star for- 
mation rate is centrally concentrated and can locally exceed 
100MQpc~ 2 Gyr~ 1 . At later evolutionary phases the gradi- 
ent flattens and star formation occurs over a large range in 
radius. In Fig.|7|we compare the time evolution of the infall 
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Figure 6. Star formation rate (upper panel) and infall rate 
(lower) versus radius. Data are shown at 2, 4, 6, 8, 10, 12, and 
13.6 (thick line) Gyrs. 
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Figure 7. Infall rate of the model versus time at different radii. 



rate at different radii. The evolution of the total star for- 
mation and infall rate of the model galaxy is shown in the 
upper plot of Fig. [g] Both distributions are very broad indi- 
cating almost constant rates over a long period of time. The 
star formation rate had its maximum ss 3Gyrs ago which is 
in reasonably good agreement with star formation histories 
of galaxies of a similar mass (peaks between 2 and 3 Gyrs) 
determin ed by the analysis o f the 'fossil record' of nearby 
galaxies jHeavens et al.ll20 04D . The present day global value 



is 3.6M Q pc-^Gyr- 



The infall rate onto the galaxy over the last WGyrs 
in the range of 2 — AM^yr -1 varies only weakly and has 
a present day value of 2.2Moyr~ 1 . This is in excellent 
agreement with estimates of infall rates of 1 — 3M0j/r _1 
that would be able to compensate for the rapid molecular - 
gas depletion in the Milky Way disc (lLacev fc FaWl98Sft . 
Observationally, a strong indication for the importance of 
infall is the signific ant presence of Deuterium at the so- 
lar neighbo urhood l|Lin5kyj |2003) as well as at the galac- 
tic centre iLubowiclie^uT |2000*) ■ As Deuterium is de- 
stroyed in stars and as there is no other known source 
of deuterium in t he Galaxy it is cosmological and of ex- 
tragalactic origin llOstriker fc Tinslevlll975l ; iLubowich et alJ 
2000). Chemical evolution models that have followed the 
evolution of Deuterium find values consistent w ith ob- 
servations llTosi et al.lll998l : IChiappini et aLlEool) . A sig- 
nificant fraction of infalling low metallicity gas could be 
high velocity clouds (HVCs) that have been accreted 
onto th e Milky Way disc and supply the ISM with Deu- 
terium dSembach et alJ [2004). Inferre d HVC infall rates 
of between G^M^yr' 1 and 5M P> yr~ 1 JWakker et alJll99Et 
iBlitz et alJll999t Braun fc Thilkerl 120041) are also in good 
agreement with indirect evidence from recent chemical evo- 
lution models that imply an average (and constant or 
eve n slowly rising) acc r etion r ate of ~ 1 — 2Mp l yr ~ 1 (see 
e.g. iGiovagnoli fc Tosil lll995l) : ICasuso fc Beckmanl d2004l) 
and references therein). This is particularly interesting 
as other successful models for th e Milky Way use ex- 
ponentially decreasing in fall rates llChiappini et alJ Il997l : 
iBoissier fc Prantzoslll999l) . 

Another independent estimate of the global infall 
rate based on the type II supern ovae rates of the Milky 
Way of ss 2 — 5century ~ 1 (e.g . [Cap j3eUjxo_el 1 _alJ ^93; 
I van den Bergh fc McClurd Il994l iTammann et alJ Il994l : 
bragicevich et alJll99s|r would result in a similar accretion 
rate for a Salpeter IMF JCasuso fc Beckmanl l2004h . This 
is in good agreement with the SNII rate derived from the 
model of w 2century~ 1 . 

The present day model values of the local SFR and IFR 
surface density are 6AMopc~ 2 Gyr~ 1 and 3.3MQpc~ 2 Gyr~ 1 
(lower panel in Fig.|HJ. For the last SGyrs the star formation 
rate varies very little. In general the local SFR is in good 
agreement with predictions from cosmologically motivated 
models for disc formation following in detail the mass aggre- 
gation histo ry of dark matter halos similar to the one of the 
Milky Way feirmani fc Avila-Reesell2000l: iHernandez et alJ 

The absolute value of the star formation rate in the solar 
neighbourhood is very difficult to deter mine observationally 
and r anges from 2 — 10M ( ? ) pc~ 2 Gyr~ 1 iGuesten fc Mezgerl 
1982). For the model star formation at the solar ra- 
dius s tarted w WGyrs ago. Observationally, ICarraro et alJ 
( 1999) have determined a similar lower limit of the age of 
the galactic disc o f 9 — l OGyrs from the ages of open clus- 
ters. |Binne^^^L| (12000,) used Hipparchos observations of 
main-sequence and sub-giant stars determine the age of the 
solar neighbourhood to be 11.2±0.75 Gyr with a lower limit 
of 9Gyr. 

The ratio of the star formation rate to the mean 
star formation rate is easier to measure than the abso- 
lute star formation rate. According to the analysis by 
iBinnev et alJ |200pj) the star formation rate at the solar 
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neig hbourhood has to hav e been very near l y cons tant (see 
also lHavwood et"allll997tl . iBertelli fc Nasi (1200 ll) found a 
star formation rate broadly increasing from past to the 
present time by comparing synthetic colour-magnitude di- 
agrams to observation of field stars from the Hipparchos 
catalo gue. Based on chromospheric activity in late type 
dwarfs lBarrvl 1^988) has found first evidence for a secular 
increase in the local star formation rate. Recent observa- 
tions of the solar neighbourhood in deed favour a fluctuat- 
ing history of the star formation rate with a slow secular 
increase. Observations of t he star formation rate history by 
iRocha- Pinto et aD I^OOOaT ) based on the chromospheric age 
distribution of nearby stars are shown in Fig. [5] For com- 
parison to our model we discarded all the stars that have 
observed ages older than the age of our model universe and 
rescaled the mean SFR accordingly (Helio Rocha-Pinto, pri- 
vate communication) . Imagine that the oldest stars also have 
the largest errors in the age determination. The observed 
fluctuations on small time scales (0.5 Gyrs) over the recent 
3Gyrs could have been trigge r ed by the passage o f spira l 
waves jHernandez et al.ll2000D . iRocha-Pinto et al.1 §000b) 
suggested close encounters with the the Magellanic clouds 
to account for the intermittent nature of the star formation 
rate and on scales of 1-3 Gyrs. The dashed line in Fig. |U] 
shows the predicted star formation rate history at the solar 
radius. To mimic the observational error we convolved the 
model data with an error of O.ldex in age (solid line). It is 
obvious that our model can not reproduce the intermittent 
nature of the SFR as none of the possible physical origins 
are included. Globally, however, our model follows the ob- 
served trend for a secular increase in the SFR history very 
well. 

In Fig. I1UI we show the mean age of the stars versus 
radius at the present day. The mean stellar age of the disc is 
5AGyrs in go od agreement with the values of the thin disc 
of 5 - IGyrs llRobin et al.ll2003|) and the mean age of the 
stars at the solar neighbourhood is about 4Gyrs. 

4.3 Metallicity evolution 

Every successful model for the formation of the Milky 
Way has to be able to reproduce the metallicity distri- 
bution of stars at the solar neighbourhood. G and K 
dwarfs, with lifetimes comparable or longer than the age 
of the galactic disc, are good tracers of the star forma- 
tion history and metal enrichment. It is now well es- 
tablished that there is a local dearth of low metallic- 
ity G- and, especially, lower mass K-stars, e.g. there are 
less low metallicity disc stars than predi cted by simple 
close d box models for chemical evolution ([van den Berghl 
119621: iPaeel fc Patchetd Il975l : iTinslevI ll980D . In particu- 
lar the distribution shows only a few stars at [Fe/H] = 
—0.4, rises stee ply to a peak at \ Fe/H ] ~ —0.1 and falls 
off t h ereafter IPagel fc Patchetl Il975l: iBeckman fe Paeell 
1 19891: ISommer-Larsenl Il99lt IRocha-Pinto fc Maciell 11996 
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proposed a revised version of the local metallicity distribu- 
tion with the peak shifted to solar metallicity. Infall of low 
metallicity gas on long time scales has long been su ggested 
as a possible solution fo r the 'G-dwarf problem' ([Larson! 
119721: iLvnden-Belll Il975t iPaeel fc Tautvaisiend Il995l) . As 
mentioned before the infall hypothesis is particularly inter- 
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Figure 8. Upper panel: Evolution of the global star formation 
rate (solid) and the global infall rate (dotted). Lower panel Lo- 
cal star formation rate (solid) and infall rate at the solar radius 
(dotted). The e rror bar indicates the obs erved range for the star 
formation rate iOuesten fc Mezge Jllflsl) . 
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Figure 9. Ratio of star formation rate over mean star formation 
rate. The dashed line shows the raw model data, the solid line is 
the model data assuming an error in the age determination of 0.1 
dex. The observed values with error ba r s are s hown by the black 
symbols. Data fromlRocha-Pint o et all l2000al) have been kindly 
provided by Helio Rocha-Pinto. 
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Figure 10. Mean age of the present day stellar disc versus radius. 
The mean age of the disc is 5AGyrs, the mean age at the solar 
radius is AGyrs. 



esting as it offers a natural explanation for the existence of 
deuterium in the galactic disc which is destroyed by star 
formation iOstriker fc Tinslevll975tlWakker et alJll99Sh . It 
has been suggested that the infall of high velocity clouds 
might continuously supply the galactic disc with metal poor 
gas. Recent observations have shown that the relatively high 
metallicities of HVC of about 0.1 s olar rule out primor- 
diallity as well as a galactic origin JSembach et al.l 12004 
IWakker et afll2003ft . 

To follow the metallicity evolution of our model galaxy 
we scaled the effective yield in Eqn. 1361 to 0.1 dex below 
the solar value at the solar radius w 4.5Gyrs ago. Implic- 
itly we assume that either the Sun has formed closer to 
the centre where the metallicity was slightly higher and 
moved outwar d thereafter probably due to orbital diffusion 
iWielenlll977l) or that the Sun had a slightly higher than 
typical metallicity representative of the ISM at the time 
of its formation. A further uncertaint y is the exact chem- 
ical composition of the Sun. Recently I Asplund et alJ <l2004l) 
have derived the solar oxygen abundance using a 3D time- 
dependent hydrodynamical model of the solar atmosphere. 
Their value of log £fo) = 8.66 ± 0.05 for the oxygen abun- 
dance and Z — 0.0126 for the solar metal mass fraction 
is significantly lower than the most commonly used value 
of logep = 8.93 and, resp ectively, Z — 0.0194 derived by 
lAnders fc Grevessd il989lt . The lower value for the solar 
metallicity is now i n better agreement with the m etallic- 



ity of the local ISM llMever et aljll99l lAndre et alJ 
nearby B stars (ICunha fcLarriber3 [^9j nKilian^t^l 



2003), 



|l994h 



and nearby young G and F dwarfs Tsofia fc Meverl kOQl). 

We used the solar metallicity value of Zq — 0.0126 
to scale our model to 0.1 dex below the solar metallicity 
Z — Zq x 10~ 01 = 0.01 at solar radius 4.5Gyrs ago, result- 
ing in an effective yield of Y = 0.0135. The effective yield 
is significantly lower than in most previous studies as we 
used the revised lower solar metallicity value for normalisa- 
tion. For the traditional value of Z — 0.0196 we would have 
needed an effective yield of Y = 0.02 if normalised to 0.1 dex 
below solar or Y — 0.025 if normalised to solar metallicity 
4.5Gyrs ago. 

Fig. II II shows the radial metallicity distribution of the 
ISM of the model disc. The distribution after 2, 4, 6, 8, 10, 
and YlGyrs is indicated by the dotted lines. For compari- 
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Figure 11. Radial metallicity distribution of the gas after 2, 4, 6, 
8, 10, 12 and 13.6 Gyrs (thick solid line). The two dotted straight 
lines indicate a slope of — 0.07dezfcpc _1 and — 0.04dexfcpc -1 . The 
metallicity gradient has been steeper in the past. The present day 
metallicity gradient of the stars is shown by the thick dot-dashed 
line. 
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Figure 12. Time evolution of the mass weighted mean metallicity 
of the stars (dashed line) and the gas (solid line). After z = 1.5 
the stellar metallicity does not show a strong evolution. 



son we also show the present day metallicty gradient of the 
stars. At present time the metallicity gradient of the model 
is about d\og(Z)/dr — 0.046dexfcpc _1 at the solar radius. 
The metallicity gradient has been significantly steeper in the 
past (especially in the inner parts of the disc) which is a typ- 
ical result for models using the inside-out formation scenario 
and has been proposed by several aut hors (Molla et al. 1997; 
Portinari fc Chiosil Il999l iHou et alJ [20001: IChiappini et alJ 
2001). Observationally, the existence of a gradient is well 
established. There is, ho wever, considerable sca tter in the 
published values (see e.g. IChiappini et alJ (|2001) for a dis- 
cussion). Our model value is on the lower end of the ob- 
served abundance gradients of light elements, X, that are 
in the r ange of -0.04 < d\oe ( X / H) / dr < -Omdexk pc' 1 
(see e.g lRolleston et ail i2000| ) : | C hiappin i et al] i 200ll) and 



references therein). Recently lEstebanetHrTf 
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Figure 13. Age-metallicity distribution of the stars at the solar 
radius for infalling gas with metallicities Zip = 1 X 10 -4 (dashed), 
Zi F = O.IZq (dotted) and Z w = Z (O.3 X 10~° 12z ) (solid). The 
effective yield has been normalised to log(Z J Zq) = —0.1 A.5Gyrs 
ago , indicated by the asterisk. Observations shown for comparison 
ajejF^ocha-PjntQ etalJ i200Cl> - diamonds with error bars - and 
lNord^tr5ni^t"^ri200^) - solid dots with typical errors shown in 
the lower part of the figure. Note that we did not fold the model 
data with the observational error in this case. 



termined the Oxygen gradient to d \og(X/H) jdr < — 0.044± 
O.OlOdeirfcpc -1 In Fig. I12l we show the time evolution of the 
mass weighted mean metallicity of the stars adn the gas 
at all radii where the total surface density is larger than 
IMq/pc 2 . As we have continuous infall of low metallicity 
gas and the high metallicity gas is depleted at the center 
for most of the time the stellar metallicity is larger t han the 
metallicity of the gas (see e.g. lAkerman et al.ll20oj) . After 
z — 1 the metallicity of the stars does not change signifi- 
cantly. The total mass of the stars, however, does increase 
by more than a factor of two (see Fig. |3J . 

The evolution of the metallicity of the ISM at the so- 
lar radius, which is equivalent to the age-metallicity distri- 
bution of the stars in the model, is shown in Fig. 1131 We 
have computed the evolution for three different assumptions 
about the metallicity of the infalling gas Zw. (a) constant 
low metallicity not pre-enriched Zip = 1 x 10~ 4 , (b) con- 
stant low metallicity enriched Zw = 0.1 Zp, which is c onsis- 
tent with HVC and Lya absorbers jBlitz et al.lll99gTl . and 
(c) a time dependent metallicity taken from cosmological 
simulations, Zip = 0.3 Zp)10~ 012z , where z is the redshift 
JCen fc Ostrikerlll999h . assuming a present day metallic- 
ity of the infalling gas of O.3Z . In general the late evo- 
lution is i n good agreement with th e age-metallicity distri- 
bution of iRocha-Pinto et alJ <l2000l) . For ages larger than 
8Gyrs the models with pre-enriched infall are more con- 
sistent wi th observations. The ag e-metallicity relation pub- 
lished bv INordstrom et al.l <|2004I) indicates an almost con- 
stant mean metallicity at all ages which might be a challenge 
to the model, especially for very old stars (for ages larger 
than 2.5Gyrs dots in Fig. 1131 . We note, however, that we 
have only plotted the raw distribution and did not take er- 
rors in the time nor in the metallicity determination into 
account. 

In the upper panel of Fig. 1141 w e show the observe d 
metallicity distr ibution of IRocha-Pinto fc Maciell il996l) . 
Ij0rgenser1 teOOOl) . and INordstrom et alJ (I2004T) . In general 



the distributions peak at values of —0.2 to —0.1 and shows 
few or no stars below [Fe/H] — — 1. The dotted thin lines 
in the lower panel of Fig |14l show the raw metallicity dis- 
tributions of the model galaxy, binned according to obser- 
vations, for the three different metallicities of the infalling 
gas. The distributions are very similar and there is a very 
good agreement between the theoreti cal and the most re - 
cently observed distributions (see also lKaratas et afll2005f) . 
Therefore pre-enrichment does not influence the local metal- 
licity distribution in our model. The constant gas infall over 
almost WGyrs determines the shape of the distribution. 
For a better comparison to the observed distributions we 
folded the average model data with a Gaussian with a dis- 
persion of 0.15 dex, which is typical for observations by 
IRocha-Pinto fc Maciel| <|l996l). an erro r of 0. 1 dex would be 
more typical for the lNoMstj5m et alJ (120041) data. The re- 
sulting distributions are shown by the solid and dashed lines 
in the lower panel of Fig. 1141 Independent of the assumed 
error there is a very good agreemen t with observations of 
the distribution of \FelH\. [rlavwoodl (1200111 has proposed a 
revision of the metallicity distribution resulting in a shift of 
the peak value to [Fe/H] — (shown as a mass fraction in 
the lower panel of Fig. 1141 . 

4.4 Photometric evolution 

We use the metal dependent models for the spect r al evo - 
lution of stellar populations of iBruzual fc Charlod {2003), 
kindly provided by the authors, assuming a Salpeter IMF to 
compute the photometric properties of our model galaxy. 

Fig |15l shows the time evolution of the absolute mag- 
nitudes, the total luminosities and the colours of the 
model galaxy at different wavelengths. The total luminosi- 
ties were calculated assuming absolute magnitudes for the 
Sun of \K, Hj R, V, B U] = [3.28, 3.32, 4.42, 4.83, 5.48, 5.61] 
jBinnev fc Merrifieldlll998T) . The absolute U and B magni- 
tudes stay relatively constant after QGyrs reflecting the al- 
most constant global star formation rate (Fig. [HJ . At longer 
wavelengths the luminosity is slowly increasing due to the 
mass assembly in the disc. The H-band magnitude of Mh = 
-23.0 (L H = 3.5 x 1O 1O L , H ) is in excellent agreement with 
the value expected from the H-Band Tully-Fisher relation 
of M H = -22.84 assuming v c = 210fcms _1 JPierce fc Tullvl 
1992). A direct comparison to the luminosity of the Milky 
Way is more complicated as the photometrical determina- 
tion of the total lu minosity of the Galaxy is v ery difficult. 
As a possible range iBinnev fc Merrifieldl Il998l) refer to ob - 
served values of L K = 4.9 x 10 10 L Q , K llKent et al.lll99lfl . 
which is close to the mode l value of Lk =4.1 x 10 10 Lq i k, 
and L L = 2.2 x 1O 1O L , L JFreudenreichlll998t) . The values 
hav e been corrected using the solar magnitudes given above 
(see lBinnev fc Merrineldlll998l) . The absolute V-band lumi- 
nosity of the model is Lv = 1.9 X lO lo L ,v which is within 
the observed range o f 1.4 x 10 10 < L v < 2.1 x 1O 1O L , V by 
Ivan der Kruhl (:1986) and slightly higher than the observed 
unobscured value of L v = 1.25 x 10 10 L Q , V (Afv = -20.4) 
(Bah call fc Soneiral fi980'). Note that uncertain dust correc- 
tions are important in fixing the observed value. The sur- 
face brightness profiles at different wavelengths are shown 
in the upper panel of Fig. 1161 Two observed values in 
the K- and B-band at th e solar radius are overplotted 
jBinnev fc Merrifieldl Il998l) . The middle panel shows the 
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Figure 14. Upper panel: O bserved m etallicity d i stribu tions of 
Rocha-Pinto & Macic 1 <1996ft ( dashed). ITorgensenl feOOOl) (solid) 
and lNoMstr5rri^T^^ri200^r 7dotted) . Lower panel: The thin dot- 
ted lines show the (almost identical) binned mean raw metallicity 
distribution of model galaxy for the models with the three dif- 
ferent infall metallicities (See Fig. 1131 . The solid line shows the 
average binned model distribution assuming an error of 0.15 dex 
in the metallicity determination. 

corresponding luminosity density profiles using the solar 
magnitudes given above. All profiles are exponential in the 
outer parts and flatten out towards the centre. The present 
day colour profiles are shown in the lower panel of Fig. 1161 
The galaxy get redder at the centre which is a combined ef- 
fect of increasing metallicities (Fig. llH and increasing stellar 
ages towards the centre (see Fig. I1UI . The B — V colour at 
the solar radius is B — V = 0.4 as opposed to an obs erved 
value of B - V = 0.6 - 0.8 feoissier fc Prantzoslll999l) . The 
B — K colour of B — K = 3.0 is more similar to observa- 
tions of B - K = 3.13 iBinnev fc Merrifieldll998T) . Over all 
the model galaxy seems to be slightly bluer than the Milky 
Way. This might be due to uncertain dust corrections. Ad- 
ditionally, the blue luminosity is very sensitive to the very 
recent (lGyr) star formation history. As we have seen from 
observations (Fig.|5J the recent star formation rate has been 
below average. This effect alone can account for the differ- 
ence in the blue luminosity of the model and the Milky Way. 



Figure 15. Upper panel: Absolute magnitude of the model 
galaxy versus time. Middle panel: Evolution of the total luminosi - 
ties with time. Observed va lues are from: fa) [K ent et al.| l|l99ll) 
. (b) iPierce fe Tullyl <1992l) . (c) iBahcall fc Soneiral Jl980l) . (d) 
Ivan der Krufd <1986ft . 



To calculate the scale lengths in the different bands we 
have fitted an exponential in the range of 1.5 < < 3.0 
scale lengths of the total surface density distribution. The 
scale lengths increase to shorter wavelengths (Fig. 1171 . This 
effect is weaker at earlier times. At present the B-band scale 
length Vd.B = 3.7kpc is a factor of « 1.2 larger than the 
K-band scale length rd,K = 3.0kpc . This trend is observed 
and is in good agreemen t with other Milky Way type spiral 
galaxies. Ide Jona Jl996T) investigated a sample of 86 nearly 
face-on spiral galaxies and concluded that spiral galaxies of 
the same type as the Milky Way have disc scale lengths that 
are a factor of 1.25 ± 0.25 larger in the B-band than in the 
K-band. This effect reflects the over all inside out formation 
pro cess for disc galaxies and ha s been found in similar stud- 
ies feoissier fc Prantzosl il999l) have found rd,B/^d.K = 1.5 
which is also in good agreement). 

In Fig llSl we show the time evolution of the K-, V-, and 
B-band surface brightness. The discs had steeper profiles 
and higher values for the central surface brightness in the 
past. This is in good agreement with recent observations of 
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Figure 16. Photometric properties versus radius. Observed 

values are = 20.62 (star) and /^b = 23.75 (diamond) 

Binncv & Mcrrificld 1998) The observed colors are B — K = 3.13 
Binncv & Mcrrificld 1998) indicated by the star in the lowest 

panel and 0.6 < B — V < 0.8 (see iBoissier fc Prantzoslll999l) 

indicated by the error bar. 

disc evolution wit h redshift that favou r the inside-out disc 
formation process ijBarden et al]l2005l) . 

Given the luminosity distributions we can compute the 
variation of the mass-to-light ratio at different wavelengths 
versus radius (Fig lliJl . In the K-band the value at the so- 
lar neighbourhood is close to 0.8 if we include the mass of 
the gas and the stars. Assuming a column density of dark 
matter within 0.5kpc of UMqpc 2 the K-band mass-to-light 
ratio would be around in ag reement with the observed value 
iBinnev fc Merrifieldl l)l998t) . However, the observed value is 
uncertain as the quantities used to calculate the mass-to- 
light ratio itself have already large errors. 



5 THE EFFECT OF A DIFFERENT IMF 

In this section we investigate the effect of a different IMF 
on the propertie s of our model disc galaxy. We tested the 
IChabrierl J2003I) IMF which provides a better fit to low 
mass stars and brown dwarfs in the Milky Way using the 
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Figure 17. Exponential scale lengths of the surface brightness 
distribution measured at the different wavelengths versus time. 
At shorter wavelengths the disc has a larger scale length. 



corres ponding evolutionary models of iBruzual fc Chariot! 
(2003). For the Chabrier IMF the fraction of gas returned 
to th e ISM by evolved stars is higher than for the Salpeter 
IMF dBruzual fc CharlotjIiooS) . Therefore we have changed 
the parameters in Erin. 1531 to R, = 0.44, <5* = 1.5, and 
tq = 1 x 10 8 to fit the corresponding return fractions of 
IBruzual fc Charlod (120031) . We have run the model using the 
Chabrier IMF with exactly the same input parameters ex- 
cept the effective yield. It has been reduced to Y = 0.013 
to guarantee a metallicity of 0.1 dex below the solar value 
at the solar radius ~ 4.5Gyrs ago. We have summarized 
the comparison in Figs. 1201 and 1211 The total gas mass and 
the local gas surface density are about 20% higher due to 
the higher gas return fraction of the Chabrier IMF (two up- 
per panels in Fig. 1201 . The stellar masses are reduced by 
the corresponding amount. Similarly the present day local 
star formation rate is about 20% higher (third panel in Fig. 
12111 and the mean age of the stars at the solar neighbor- 
hood is reduced to 3.8Gyrs. The effect on the luminosi- 
ties is stronger as the Chabrier IMF in general leads to 
1.4 -1 .8 times smaller mass-to-light ratios than the Salpeter 
IMF dChabrierll20"03l: IBruzual fc CharlotJl2003l) . e.g. the ri- 
band luminosity increases from Lk = 4.1 x 10 10 Lk© to 
L K = 7.2 x 10 10 Lkq (fourth panel in Fig.^OJ- The lowest 
panel in Fig. I20l shows the time evolution of the local colors 
that are only weakly effected by th e different IMF as already 
mentioned bv IBruzual fc Charlotj ll2003l) . 

In the upper panel of Fig. 1211 we show a comparison 
between the radial metallicity distribution of the ISM for a 
Salpeter and a Chabrier IMF. The Chabrier IMF leads to a 
higher metallicity at the central region an the present day 
gradient at the solar radius is increased to — 0.056de:rfcpc _1 . 
The Chabrier IMF results in a larger fraction of stars with 
super solar metallicity at the solar radius but the shape of 
the metallicity distribution at the solar radius remains ba- 
sically unchanged (lower panel in Fig. I2H . 

Taking the observational errors into account the proper- 
ties of the model using a Chabrier IMF are as well in agree- 
ment with observations, especially the metallicity gradient 
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Figure 18. Evolution of the K-, V-, and B-band rest frame sur- 
face brightness profiles with time. The dashed line from the left 
to the right show the surface brightness profiles from early times 
to now. The profiles are plotted every lAGyrs. The present day 
profiles is the solid line (see Fig. 1161 . 
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Figure 19. Present day disc mass-to-light ratio (stellar + gaseous 
mass, the exact dark matter contribution is unknown) versus ra- 
dius at different bands. Dark matter is not included here as the 
exact distribution is unknown. Assuming a column density of dark 
matter within 0.5kpc of IIMqjjc 2 at the solar radius the K- and 
B-band mass-to-light ratios would be 1 and 1.6, respectively. 
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Figure 20. Comparison of the time evolution of already intro- 
duced characteristic model parameters for a Salpeter IMF (thin 
lines) and a Chabrier IMF (thick lines). 
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Figure 21. Upper panel: Present day radial metallicity distribu- 
tion of the gas for the Salpeter IMF (thin dotted line) and the 
Chabrier IMF (thick solid line). The two dotted straight lines 
indicate a slope of — 0.07dexkpc~ 1 and —O.OAdexkpc -1 . Lower 
panel: Local metallicity distribution of the stars for the Salpeter 
IMF (thin dotted line) and the Chabrier IMF (thick solid line). 



and the local metallicity distribution. The local value of the 
gas surface density might be at the upper limit of observa- 
tions (the gas distribution at the solar radius is very difficult 
to measure) but could be reduced by assuming a higher star 
formation efficiency. Therefore we can not conclude that ei- 
ther the Salpeter IMF or the Chabrier IMF leads to model 
properties that are over all in better agreement with obser- 
vations. 



THE EFFECT OF STAR FORMATION 
THRESHOLDS 



iMartin fc Kennicuttl (1200 ll) have shown that disc galaxies 
show prominent breaks in their outer Ha profiles giving ev- 
idence for the existence of a critical surface density beyond 
which the star formation is strongly suppressed. A possible 
consequence of a star formation threshold density could be 
the obser ved truncation or outer breaks in the pro file of stel- 
lar discs Ivan der Kruidl979l lde" Griis et al However, 



other possible explanations, e.g. trunca tion by the maximum 
angular momentum of the protogalaxy jvan der Kruitll987l) , 
have been given as well. As we have seen in Section |5] there 
is observational evidence that the stellar disc of the Milky 
Way is truncated at a radius of « 15fcpc. In this section 
we investigate whether this truncation can be caused by a 
threshold gas surface density for star formation. 

Based on observational constraints any star formation 
law for disc galaxies will break down at gas surface densities 
of 



Ecrft < 5 - 10M Q pc 



(37) 



where the star f ormation is in general strongly suppressed 
iKennicutd Il998l) . However, in individual cases the ac- 
tual threshold value may vary by an order of magnitude 
jMartin fc KennicuttihoOll) and the simple prescription is 
not valid. It has been argued for a long time that gravita- 
tional instabilit ies could determine the threshold density for 
star formation jToomrelll96l lQuirklll972l ; lKennicutd[l989l) 
by coriolis forces 



E C rit = QQ ■ 



7tG 



(38) 



where <T ga s is the velocity dispersion of the gas, k is the 
epicyclic frequency, 



- M = 2(- + -- 



(39) 



and qq = 1 IQ defines an effect ive Toomre stability parame- 
ter. [^^rth^^Kgmii^ii^] l)200lT) favour a value of «q = 0.69 
corresponding to Q = 1.45 for their sample. Numerical stud- 
ies of self-gravitating g a s discs indicate v alues up to Q — 1.6 
<Kim fc Ostrikerll200lt iKim et al]l2003l) Alternatively, the 
disc can be stabilised by shear forces 



S C rit = «A- 



,-1 



7tG 



(40) 



where A is the Oort constant iHunter et al.lll99ci) . For con- 
stant circular velocities v c like in the outer parts of disc 
galaxies both criteria reduce to 



Ecrit = Q!Q 



V2o 



ivGr 



(41) 



2-kGt 

and oq = 0!a/(2\/2)- Both criteria differ in the innermost 
part of the disc if the circular velocity is not constant. 

The above criteria only consider the theoretical case of 
a purely gaseous disc. In real galaxies the s tellar compo- 
nent c an have a destabilising effect on the gas. lWang fc Silkl 
lll994l) derived an effective stability parameter a e a taking 
the stellar surface density E* and velocity dispersion a, into 
account 



Ecrit 



(JK 

«cff —p, , where a cff 

7I~G 



E*<To 



(42) 



They have successfully applied their criterio n to the M ilky 
Way, which has recently been confirmed by iBoissier et alJ 
J2005D . 

We have tested the influence of the threshold prescrip- 
tions given by Eqns. 1371 1411 and 1421 on the evolution of 
the outer parts of our model disc. For the constant thresh- 
old case (Eqn. 1371 we assumed a critical gas density of 
E C rit = 7Mqpc~ 2 , for the Toomre based criterion (Eqn. HT1 
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we used the observed value of «q = 0.69 and a velocity dis- 
persion of the gas of <7 gas = 6fcms _1 . For the third version 
we assumed a stellar velocity dispersion based on Q =1.6 
to compute a e g. We mimicked the fluctuating nature of the 
gas surface density at a given radius by randomly drawing a 
surface density from a Gaussian with a mean of the model 
surface density and a dispersion of 10% that is compared to 
the threshold value. This procedure guarantees that some 
star formation can occur even in regions where the mean 
surface density is below threshold. As soon as the gas den- 
sity drops below the critical density we let the star formation 
rate cease immediately 



1000.0 



£sfr(M) =0 for E gas < E c 



(43) 



We focus on effects on the outer disc at radii larger 
than 7kpc. In the central parts a simple threshold prescrip- 
tion is less meaningful due to the limitations of some of the 
prescriptions (see above) and the uncertain influence of a 
central bulge or bar component which is not explicitly in- 
cluded in our model. The effect of the cutoff on the outer 
mass, gas and the K-band light distribution is shown in Fig. 
1221 At radii smaller than « 12kpc, the cutoff has almost no 
effect. At radii of about 12 — 13fcpc the cutoff models produce 
a sharp decline in the stellar distribution in agreement with 
observations of the Milky Way and external disc galaxies. 
The effect on the K-band light profile is shown in the low- 
est panel of Fig. 1221 All models produce a similar behaviour. 
The gas surface density is slightly increased compared to the 
model without threshold at radii larger than 12kpc and is 
qualitatively in better agreement with observations. Other 
physical mechanisms might be im portant for disc trunca- 
tions as well (see e.g. Schavc 2004). However, based on this 
model we conclude that even a simple recipe for star for- 
mation cutoff can be successfully implemented in analytical 
models and producing truncated present day stellar discs in 
good agreement with observations, Qualita tively our results 
are in agree ment with previous models bv IChiappini et ail 
il997l l200lf) who used fixed threshold value. At the solar 
neighbourhood a cutoff leads to a delayed onset of star for- 
mation (Fig. I23(l . For all models the star formation sets in 
after ss 5 — 6Gyrs. 
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7 SUMMARY & CONCLUSION 

We have presented a simple model for the evolution of the 
Milky Way in a cosmological context. We start with the 
current observed gross properties of the Galaxy. Assuming 
simple scaling relations derived from the standard properties 
of a concordance ACDM cosmological model that determine 
the time evolution of the infall rate onto the galaxy and also 
the time dependence of the rotation curve we have our back- 
bound state. We supplement this with a standard star for- 
mation prescription, and chemical evolution models assum- 
ing a Salpeter IMF and are able to reproduce a large number 
of observed global and local properties of the Milky Way to 
a similar accuracy as previous successful models for chem- 
ical evolution and chemo-spectrophotometr ic models (e.g. 
IChiappini et~allll997t IPrantzos fc SiUdll998l and references 
therin). However, in contrast to most previous models the 
evolution of the infall rate itself with time and radius is not 
a free parameter in the model presented here. 



Figure 22. The stellar surface density (upper panel), the gas 
surface density (middle panel) and the K-band light distribution 
(lower panel) without star formation cutoff (solid line) and with 
cutoff for three different prescriptions for the critical surface den- 
sity: constant at 7Mqpc~ 2 (dotted), based on the instability of 
the gas disc (dashed) and based on the instability of the combined 
stellar and gaseous disc (dashed-dotted) . 

Globally, the gas infall rate and the star formation rate 
are almost constant with time at w 2 — 4M0j/r _1 during 
the evolution of the disc. The total mass of the present 
day model galaxy in stars and in gas is 2.7 x 10 10 A/q and 
9.5 x 10 9 Mq, respectively, and the disc has an absolute ri- 
band magnitude of —23.2. The present day scale length at 
longer wavelengths is shorter than at shorter wavelengths. 
The mean age of long lived stars at the solar neighbour- 
hood is about 4Gyrs. Presently, the local surface density 
of the stars and gas are 35 and 15Mqpc -2 , respectively. 
We have been able to reproduce a narrow metallicity dis- 
tribution of the stars at the solar radius with a peak at 



A simple model for the evolution of disc galaxies: The Milky Way 19 




Figure 23. Star formation history at the solar neighbourhood for 
different models of cutoff. In general, star formation cutoff leads 
to a later onset of star formation and a slightly younger age of 
the solar neighbourhood, (see Fig. I§1 



[Z/Zq] = —0.1 and a radial gradient of — 0.046dea;fcpc _1 
which has been significantly steeper in the past. We have 
tested different threshold density prescriptions for star for- 
mation that all lead to a truncation of the stellar disc at 
about Ylkpc. Changing the IMF from Salpeter to Chabrier 
leads to an increased gas content and higher luminosities of 
the model galaxy. The present day metallicity gradient is 
steeper and the local metallicity distribution is only weakly 
affected. Based on our results we can not decide whether a 
Chabrier IMF or a Salpeter IMF leads to a better agreement 
with observations. In total we have been able to produce well 
the local age and metallicity distributions and the star for- 
mation history which provide the local fossil record of our 
Galaxies' evolution. 

The model described here is, of course, over simplified. 
Several physical processes that are likely to play a role dur- 
ing the formation of disc galaxies have been excluded. In 
the following we mention only a few processes that might 
have an influence on the over all results: In the hierarchical 
framework, disc galaxies do not only grow continuously by 
smooth accretion of gas but also in discrete steps by mi- 
nor mergers and accretion of small satellites. Simulations 
of disc galaxy formation have shown that not only gas is 
building up the disc but a small amount (up to 10%) of old 
stars that have formed within small satell ites are continu- 
ously added to t he disc by accretion iAbadi et al.l [2003b; 
iMeza et al.ll20o3) . A process that is presently observed in 
our Milky Way for the Sagg itarius dwarf galaxy and other 
small satellites (e.g . |jl3aJaL^L-ajJ|l997|: iNewberg et alJl2002t 
iMartin et alJ 12004 iHelrni et alJl200rJ) . Accretion of satel- 
lites might add very old stars to the solar neighbourhood 
thereby increasing its apparent age. The smooth accretion 
as assumed in the model can only be seen as an approximate 
representation of the past integrated accretion rate. In real 
disc galaxies minor mergers, satellite encounters and the pas- 
sage of spiral waves might locally enhance the star formation 
rate and lead to the intermittent star for mation history that 
is observed in the local neighbourhood iRocha-Pinto et all 



l2000al ; Iflernandez et aTll200Cl) . This oscillatory behaviour 
can not be reproduced by the model. In particular, we have 
not been able to exactly fit the observed colours at the so- 
lar neighbourhood. We have tested oscillatory star forma- 
tion rates with a local depression at the solar neighbour- 
hood about 1 Gyr ago and found a good agreement with 
the data. The present day local and global colour of discs 
can not strongly constrain any simple and generally valid 
model on disc galaxy evolution. In particular, the influence 
of an oscillatory star formation history on the metallicity 
evolution is very small as the metallicity is basically deter- 
mined by the past int egrated star formation history (see also 
iHirashita et alJl200lTl . 

In this paper we have focused on the evolution of the 
disc component only and did not include explicitly the effect 
of gas infall from a stellar spheroidal component (halo) of the 
Milky Way. This process, that has been i nvestigated by sev- 
eral authors, has been proposed (see e.g . lOstriker fc Thuara 
ll975tlRandll99ll:IWvse fc Gilmorelll992h as a possible solu- 
tion to the G-dwarf problem as it would provide pre-enriched 
gas from halo stars. Implicitly, we have investigated this ef- 
fect by allowing the metallicity of the infalling (to the disc) 
gas to have the metallicity expected from gas pre-processed 
through halo stars and the gas expelled form nearby dwarf 
systems. The present metallicity distribution at the solar 
neightborhood is determined by the almost constant gas in- 
fall rate over the last WGyrs is only very weakly effected by 
pre-enrichment. We will focus on the role of bulge compo- 
nents in the framework of our model in a subsequent paper. 

The influence of radial gas flows possibly triggered by 
tidal interactions, infalling low angular momentum gas or 
viscosity in the g as disc have been investig ated by several 
authors (see e.g. IPortinari fc Chiosl J2000) and references 
therein). If important inflowing gas is likely to increase the 
metallicity gradient and steepen the density profile. On the 
other hand orbital diffusion of stars might flatten the gradi- 
ent again. In addition, we assumed instantaneous mixing of 
recycled material within one ring and avoided complicated 
feedback and mixing processes that might be more or less 
importa nt in the early o r late phases of the disc evolution 
(see e.g. lEfstathioulkoOOl) 

All the above processes might play a role during disc 
galaxy evolution but were not implicitly included in the 
model. Still the presented results look very promising. We 
might conclude that to first order smooth gas infall followed 
by local star formation as a function of the local gas density 
is the dominant process shaping the Milky Way and other 
disc galaxies. 

In this paper we have taken an archaeological approach 
to galaxy evolution - to check if the computed history is 
consistent with the local fossil record. Alternatively, one can 
compare the evolution directly with the past history as seen 
in for instance in the HDF or HUDF or redshift surveys 
like GEMs. We will pertain to this test in a subsequent pa- 
per but simply note that the evidence is consistent with 
an inside-out formation scenario of discs, e.g. smaller aver- 
age discs sizes at earlier times. Recent high redshift obser- 
vations strongly support this scenario |Ravindjariath_e^ - ^l 



vations strongly support tins scenario Kavmdranath et al 
2004l:lFerguson et all2004l:lBarden et all2005HTruiillo et al 



2005) 
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